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The triple-a reaction rate is re-evaluated by directly solving the three-body Schrodinger 
equation. The resonant and nonresonant processes are treated on the same footing using 
the continuum-discretized coupled-channels method for three-body scattering. Accurate de- 
scription of the a-a nonresonant states significantly quenches the Coulomb barrier between 
the two-a's and the third a particle. Consequently, the a-a nonresonant continuum states 
below the resonance at 92.04 keV, i.e., the ground state of 8 Be, give markedly larger contri- 
bution at low temperatures than in foregoing studies. We show that Nomoto's method for 
three-body nonresonant capture processes, which is adopted in the NACRE compilation and 
many other studies, is a crude approximation of the accurate quantum three-body model 
calculation. We find about 20 orders-of- magnitude enhancement of the triple-a reaction rate 
around 10 7 K compared to the rate of NACRE. 

After the Big Bang, all elements surrounding us have been created in the cosmos. 
Among them, 12 C is one of the most important nuclei because it is an essential 
element of life. In this sense, to understand the origin of 12 C would be equivalent to 
know "where we come from. " In view of nuclear physics, it is well known that the 
second + state of 12 C at 7.65 MeV above its ground state was predicted by Hoyle in 
1953 to explain the abundance of 12 C, or, existing of human beings; this challenging 
prediction was soon confirmed experimentally, x ) and the 0~j~ state newly discovered 
is called the Hoyle resonance. Since this historic discovery of the Hoyle resonance, 
the a(aa, 7) 12 C reaction, i.e., the so-called triple-a reaction, has been described as 
a series of the following two reactions: 2 ) 

a + a^ 8 Be, a + 8 Be -> 12 C(2+) + 7. (1) 

In the latter, 12 C in the 2f state is formed by a 7-decay from the Hoyle resonance, 
and then decays into the ground state of 12 C, with emitting 7 again. This picture 
of the triple-a reaction, which is realized at temperatures higher than a few 10 8 K 
in helium burning stars, successfully describes the present abundance of 12 C. 

At low temperatures, however, the energy of the three a system cannot reach 
the resonance energy of 379.5 keV above the three-a threshold. A formation pro- 
cess through nonresonant three-a continuum states, i.e., the nonresonant triple-a 
process, then becomes dominant. The formation rate of 12 C at low temperatures 
is of crucial importance for studies on helium burning in accreting white dwarfs 
and neutron stars, 3 )' 4 ) and is believed to strongly affect the evolution of primordial 
stars. 5 ) In Ref. 3) Nomoto proposed a method for evaluating contribution of the 
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nonresonant triple-a process by using resonance formulae, with an energy shift for 
the Hoyle resonance, to describe the two reactions of Eq. ([1]). This method, which we 
call Nomoto's method in this Letter, was found to give a significantly larger triple- 
a reaction rate at low temperatures 3 )' 6 ) than that obtained by a naive resonance 
formula without the energy shift. Nomoto's method has been a standard method 
for describing the nonresonant triple-a reaction. It should be noted, however, that 
Nomoto's method does not explicitly describe the role of the nonresonant continuum 
states in the triple-a process. 

In the present Letter, we evaluate the nonresonant triple-a reaction rate by di- 
rectly solving the Schrodinger equation of the three-a system. The three-a scatter- 
ing wave function is obtained by the continuum-discretized coupled-channels method 
(CDCC), 7 ) which was proposed and developed by the Kyushu group more than 20 
years ago, and has been successfully applied to studies of various three-body reaction 
processes; see, e.g., Refs. 7)— 9). In CDCC, resonant and nonresonant states of the 
a-a system, as well as those of the three-a system, are treated on the same footing. 
This is one of the most important advantages of the present calculation. 

As for previous three-body calculations of the triple-a reaction, Kamimura and 
Fukushima 10 ) firstly showed, using the microscopic three-a resonating group method, 
the importance of the couplings between the a-a resonant channel and the a-a 
nonresonant channels for reproducing the properties of the Hoyle resonance and 
understanding the processes of Eq. (JTJ) . This finding was confirmed by Descouvemont 
and Baye 11 ) using the microscopic three-a generator coordinate method; they derived 
the triple-a reaction rate down to 10 7 K. In these studies, however, treatment of the 
a-a nonresonant continuum was very primitive; only a few discretized nonresonant 
states were included. More seriously, a-a nonresonant continuum states below the 
resonance at 92.04 keV, which play essential roles in the nonresonant triple-a process 
as shown below, were completely missed. Thus, it is obvious that these models 
cannot accurately describe the nonresonant triple-a process at low temperatures, 
say, T < 10 8 K. 

CDCC is well known as one of the most accurate reaction models to describe 
three-body reactions. It works even cases where Coulomb interactions play dominant 
roles and careful description of low-energy continuum states is required. In the 
present study, we use the three-a scattering wave function obtained by CDCC, and 
evaluate the triple-a reaction rate at temperatures of 10 7 K < T < 10 9 K. Emphasis 
is on the reaction rate at low temperatures for T < 10 s K, where contribution of the 
nonresonant triple-a process is dominant. We show that Nomoto's method 3 )' 6 ) used 
in the NACRE compilation 12 ) and many other studies is a crude approximation of 
the accurate quantum three-body model calculation. At T = 10 7 K, typically, our 
new result is more than 20 orders-of- magnitude as large as that of NACRE. 12 ) 

In the present calculation, we work with the three-a system shown in Fig. Q] 
based on the Jacobi coordinates (r, R). The relative energy of a\ and ai is denoted 
by ei2 and the relative energy between 03 and the center-of-masses (cm.) of a\ 
and 02 is denoted by €3; the total energy E in the three-a cm. frame is ei2 + £3. 
We obtain the triple-a reaction rate by directly solving the following three-body 
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Fig. 1. Illustration of the three-a system. 

Schrodinger equation: 

[T r + T R + v(r) + + v(R 2 ) - E] #(r, R) = 0, (2) 

where T r and Tr are the kinetic energy operators associated with r and R, re- 
spectively, and v is the interaction between two-a's consisting of the nuclear and 
Coulomb parts. 

Following the standard continuum-discretization procedure called the average 
method, 7 ^ we first discretize the continuum states of the a\-a 2 subsystem into mo- 
mentum bins. We prepare L 2 -integrable a\-a 2 wave functions iii(r) (i = l-« m ax) 
by 



i+i 



Ui{r) = I fi(k)u(k,r)dk 



\h{k)\ 2 dk 



1/2 



(3) 



where u(k, r) is the a±-a 2 scattering wave function with the relative momentum k, 
and fi(k) is a weight function set to be a constant and the Breit-Wigner function 
for nonresonant and resonant bin states, respectively. 7 ^ Normalization of u(k, r) is 
defined by J u*(k\ r)u(k, r)dr = 5(k'— k), which makes Mj(r) satisfy J u*(r)v,j(r)dr = 
5{j. The average momentum (energy) of the ith bin state of the ct\-a 2 system is 
denoted by k{ (ei2,i) hi the following. The total wave function of the three-a system 
with CDCC is given by 

& <rR) = ^ i_^w£!m 

where xi°\r) describes the relative motion between {a\-a 2 ) in Uj and with 
the relative momentum Ki that is obtained by energy conservation of the three-a 
system. The index iq represents the incident channel. Note that we consider only the 
s- waves of ui and x[ 1 °\r), since we are interested in a reaction at very low energies. 
Furthermore, we have dropped all phase factors in giving no contribution 

to the reaction probability shown below. The coupled-channel (CC) equations for 
x\ 1 °\r) (i = l-«max) are given by 

[Tr + Vu (R) — (E — ei 2)i )] X? o) (R) = V »' W ^ 
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which are solved under a usual boundary condition for xf°\R). 7 ^ 
potential Vw (R) is defined by 



Viv (R) 



Ui (r) 



v(R 1 ) + v(R 2 ) 



The coupling 



(5) 



The reaction probability {ov)- k _ E of the triple-a process due to the electric 
quadrupole (E2) transition to the 2^~ state is given by 



2(2vr) 7 (hw\ 
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where !?w is the wave function of the 2^~ state of C with M the projection of the 
total spin, is the E2 transition operator, and the photon energy fko is given 
by 2.8358 + E MeV. Note that we use the symbol [av) for the reaction probability 
following NACRE; 12 ) it is actually the E2 transition probability divided by the nor- 
malization factor l/(2-7r) 6 for the three-a scattering wave function in free space. The 
triple-a reaction rate (aaa)(T), as a function of T, is obtained by taking an average 
of (crv)^, E with the Maxwell-Boltzmann distribution for the velocity of each a: 



with 



aaa)(T)=3Nl-^- 3 J j|> M^jexp (-— j 
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Fig. 2. Qi-Q2 momentum bin states included in the present CDCC calculation. 
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where k-Q is Boltzmann's constant and N\ is Avogadro's number. Factor 3 in Eq. ([6]) 
comes from the fact that the three-a system is symmetric with respect to an exchange 
of each pair of alpha particles. 12 ) 

In numerical calculation, we discretize the fc-continuum of the a\-a,2 system from 
0.008 fin" 1 (ei2 = 0.668 keV) to 0.130 fni" 1 (e 12 = 176 keV) with the width of 0.001 
fm , which results in z max = 122. Schematic illustration of the a.\-a.2 bin states 
taken in the CDCC calculation is shown in Fig. [2l This discretization is sufficiently 
precise for the present purpose. The 86th bin corresponds to the a\-cti resonance. 
The maximum value of r is 5,000 fm with the increment Ar of 0.1 fm. 

As for the nuclear potential between two a's, we use the following two-range 
Gaussian form (with depth in unit of MeV): 

v™ c \x) = lOO.Oe"^ 1 - 00 ) 2 - 30.35e"^ 2 - 13 ) 2 , (7) 

where x is the displacement of the two a's in fm. The first repulsive part simulates 
the Pauli exclusion principle that nucleons in an a cannot occupy the nucleon s-orbit 
in the other. This potential gives a-a resonance at 92.0 keV with the width of 4.8 eV, 
which reproduce well the corresponding experimental values, i.e., 92.04 ± 0.05 keV 
and 5.57 ±0.25 eV. 13) It is known 14 ) that this type of a-a potential with a repulsive 
core is not suitable for describing the 0^ or 2f state of 12 C, in which the three-a 
particles are closely bound. However, it can successfully be applied to the 0Jj~ state, in 
which the three-a's are loosely coupled. Therefore, we consider that the simulation 
of the Pauli principle by introducing the repulsive part is justified in describing the 
three-a scattering states. In fact, it is numerically confirmed that even if we put 
v nucl (x) = in the calculation of Vat (i?) given by Eq. ([5]), the resulting reaction rate 
(aaa)(T) for T < 10 8 K changes by only about 2% at most. 

Equations (J2|) are numerically integrated up to i? max = 2, 500 fm with AR = 0.25 
fm, and xf o) (R) (i = 1- imax) ars connected to usual asymptotic form. The total 
energy E is varied from 1 keV to 500 keV with AE = 1 keV; around the Hoyle 
resonance at E = 379.5 keV, we put AE = 0.1 keV. In the evaluation of the coupling 
potentials Vu/ (R) given by Eq. ([5]), we reduce v llucl by 1.5% so that the (ai-a2)-a3 
system in the i = 86 channel, with iq = 86, forms a resonance at = 287.5 keV. 

In the calculation of tyfj , use of the potential of Eq. © is not appropriate 
as mentioned above. Instead, we adopt a sophisticated three-a wave function 15 )' 16 ) 
that was obtained on the basis of the orthogonality condition model 17 ) for describing 
the Pauli principle with the three-a particles symmetrized. In this semi- microscopic 
calculation, the a-a potential was derived by folding an effective nucleon-nucleon 
force 18 ) into the a-particle density. 

We show in Fig. [3] the discretized continuum wave functions Ui(r) of the ai-a2 
system (upper panel) and the Coulomb parts of the diagonal coupling potentials 
V i f(R) (lower panel). We use logarithmic scale for the horizontal axis in each panel. 
The solid, dashed, and dotted lines correspond to the ai-a2 continuum states with 
the average energies i\2,i of 38.2 keV (i = 53), 92.0 keV (i = 86), and 152 keV 
(i = 113), respectively. One sees that the resonant wave function has a dominant 
amplitude in the interaction region (r < 10 fm), while the nonresonant wave functions 
have appreciable amplitudes only at larger r. This clear difference in Ui{r) drastically 



G 



Letters 



affects V^{R) as shown in the lower panel. One sees that the Coulomb barrier height 
V^{R) for the 01-02 nonresonant bins is very much lower than for the resonant bin. 
Thus, 03 can easily approach to the a\-cti system when a.\-ai is in nonresonant 
states. This important feature has not been considered in preceding studies with 
Nomoto's method. 

Figure [4] shows the calculated results of {w)j t _ E ; to show the difference between 
the resonant and nonresonant results clearly, CC effects are not included here. The 
solid (dashed) line corresponds to io = 53 (86). One sees that the solid line has a 
completely different energy dependence from that of the dashed line. The nonres- 
onant reaction probability E is almost constant (in the scale of the vertical 

axis in Fig. U|) above E ~ 70 keV, and dominates the resonant one for E < 200 keV. 
Note that for iq = 53, E = 70 keV corresponds to £3 = 31.8 keV that agrees well 
with the quenched Coulomb barrier height shown by the solid line in Fig. [3] (lower 
panel). 

If we disregard unphysically the channel dependence of the coupling potentials 
Vn'(R) and take only the diagonal components, by using the following replacement: 

V W (R) ^86,86 (R)S ii/ , (8) 

we can simulate the evaluation of (av) r „ with Nomoto's method; note that the 86th 

bin corresponds to the cti-ct2 resonant state. This replacement of Va'(R) in Eq. (JU) 
makes the three-a system form a resonance when €3 = 287.5 keV, independently 
of ei2- In other words, when e±2 = 92.0 — Ae\2 keV, the three-a system has a 
resonance at E = 379.5 — Ae\2 keV. This Ae±2 is nothing but the energy shift in 
Nomoto's method. The result for iq = 53 with Eq. ([5D is shown by the dotted line in 




Fig. 3. (Color online) Discretized 01-02 continuum wave functions u;(r) (upper panel) and the di- 
agonal Coulomb potentials Vu(R) (lower panel). The solid, dashed, and dotted lines correspond 
to i = 53, 86 (resonance), and 113, respectively. 
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Fig. 4. (Color online) 
shows the result for io 
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for io = 53 (solid line) and io = 86 (dashed line). The dotted line 
53 with Eq. 1(5). which simulates Nomoto's method. 
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Fig. 5. (Color online) Triple-a reaction rate as a function of temperature. The solid line represents 
the result of CDCC. The dashed line shows the contribution of resonant capture. The result of 
CDCC simulating Nomoto's method with Eq. ([8| is shown by the dotted line. The dash-dotted 
line shows the reaction rate of NACRE. 



Fig. [H Although it has a peak at E ~ 326 keV as expected, its energy dependence is 
similar to that of the dashed line, which results in a much smaller value than the true 
nonresonant capture probability (solid line) at low energies. Thus, Nomoto's method 
is shown to be a crude approximation of the accurate three-body model calculation. 

We show by the solid line in Fig.[5]our result of the triple-a reaction rate obtained 
by CDCC. The triple-a reaction rate of the NACRE compilation 12 ) is shown by the 
dash-dotted line for comparison. Drastic enhancement of the reaction rate at low 
temperatures is found. The dashed line shows our result including only the resonant 
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capture process, i.e., via the {ol\-ol2)-o-z resonance state at E = 379.5 keV. Differ- 
ence between the solid and dashed lines clearly shows the dominant contribution of 
the nonresonant triple-a process for T < 2 x 10 8 K. As mentioned above, we can 
simulate the NACRE evaluation, which is essentially based on Nomoto's method, by 
using Eq. (J8]); the result of this calculation is shown by the dotted line. As expected, 
the dotted line reproduces well the result of NACRE. Thus, we conclude that the 
nonresonant capture process has much larger contribution than in previous evalua- 
tions, 3 ^' 6 )' 12 )' 19 ) as a result of the significant reduction of the Coulomb barrier height 
between and the nonresonant ai-c^- This barrier reduction cannot be taken into 
account if one uses Eq. (jSJ), or, equivalently, adopts Nomoto's method as mentioned 
above. Another remark is on the importance of the low-energy ol\-oli continuum 
states below the resonance at 92.04 keV, which were completely missed in the pre- 
ceding three-a model studies. 10 )' 11 - 1 It is found that these low-energy states have 
almost all contributions to the total reaction rate (solid line) for T < 4.0 x 10 K. Be- 
cause of the lack of these states, the reaction rate given in Ref. 11), e.g., 5.34 x 10 -63 
cm 6 s _1 mol~ 2 at 10 7 K, is markedly smaller than the present result at low temper- 
atures. 

Table I. The triple-a reaction rate (in cm 6 s _1 mol -2 ) obtained by CDCC, together with its ratio 
to the rate of NACRE. The number a[n] means a x 10". Rates at other temperatures are 
available. 



T 


(acta) 


ratio 


T 


(acta) 


ratio 


(10 7 K) 






(10 7 K) 






1 


1.08 [-44] 


3.7[+26] 


15 


1.52[-16] 


9.5[+01] 


1.5 


3.42[-38] 


5.4[+23] 


20 


1.92[-15] 


1.9[+0] 


2 


3.12[-34] 


5.7[+21] 


25 


4.37[-14] 


1.0[+0] 


2.5 


1.73[-31] 


1.6[+20] 


30 


4.51[-13] 


9.9[-l] 


3 


2.44[-29] 


1.7[+18] 


35 


2.29[-12] 


9.8[-l] 


4 


1.10[-25] 


2.1[+15] 


40 


7.37[-12] 


9.8[-l] 


5 


3.41[-23] 


3.3[+13] 


50 


3.41[-11] 


9.9[-l] 


6 


1.63[-21] 


1.4[+12] 


60 


8.56[-ll] 


9.9[-l] 


7 


2.56[-20] 


8.5[+10] 


70 


1.54[-10] 


9.9[-l] 


8 


2.01[-19] 


2.1 [+09] 


80 


2.26[-10] 


1.0 [+0] 


9 


9.89[-19] 


3.9[+07] 


90 


2.93[-10] 


1.0 [+0] 


10 


3.52[-18] 


1.5[+06] 


100 


3.48[-10] 


1.0 [+0] 



For more detailed comparison, we need renormalization for the CDCC result. 
We introduce an effective charge be = 0.77e to the E2 transition operator, so that 
our result of the B(E2) value, evaluated at E = 379.5 keV with the + wave function 
normalized to unity, reproduces the experimental value of 13.4 e 2 fm 4 , 20 ) just in the 
same way as in Ref. 11). Then, we renormalize our result (with 5e) to the reaction 
rate of NACRE at T = 10 9 K, where the resonant capture process through the 
Hoyle resonance is dominant; 12 ) the renormalization factor obtained is 1.54. From 
this value, one may estimate the uncertainty of the present calculation to be more 
or less 50%. 

In Table U we show the renormalized triple-a reaction rate obtained by CDCC, 
together with its ratio to the rate of NACRE, at some typical temperatures. As ex- 
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pected, for T > 2.5 x 10 K the ratio is almost unity that shows our calculation, with 
the renormalization at 10 9 K, reproduces very well the temperature dependence of 
the resonant triple-a process through the Hoyle resonance. At low temperatures, on 
the other hand, the ratio exceeds 10 20 , which will affect the helium burning in accret- 
ing white dwarfs and neutron stars. Furthermore, even at rather high temperature of 
1.5 x 10 8 K, we obtain a larger reaction rate by almost two orders-of-magnitude than 
that of NACRE. It should be noted that a broad 2;j~ resonance state at E = 1.75 MeV 
with the a-decay width of 0.56 MeV 11 ) is included in the NACRE evaluation. In 
fact, it is shown that this 2% state has significant contribution to the triple-a reac- 
tion rate for T > 2 x 10 9 K. 12 ) At these high temperatures, the nonresonant triple-a 
process that we focus on in the present study will obviously be negligible compared 
to resonant processes via the 0^, 2^~, and other possible resonance states 5 ) of 12 C. 
Very recently, it was reported 21 ' that a stellar evolution model computed with our 
new reaction rate of the triple a reaction caused inconsistency with observations of 
red giant branches. Further investigation on this implication will be very interesting 
and important. 

In summary, we have evaluated the triple-a reaction rate by directly solving the 
three-body Schrodinger equation with CDCC. We treat the resonant and nonreso- 
nant processes on the same footing. The a-a continuum states below the resonance 
at 92.04 keV are shown to play essential roles in the nonresonant triple-a process for 
T < 4.0 x 10 7 K. The key of the nonresonant capture process is that the Coulomb 
barrier between the two-a's and the third a is much quenched compared to that in 
the resonant capture. This property extremely enhances the nonresonant triple-a 
process at low temperatures, i.e., T < 10 8 K. The ratio of our triple-a reaction 
rate to that of the NACRE compilation is more than 10 26 at 10 7 K and about 100 
at 1.5 x 10 8 K. It is found that Nomoto's method for three-body nonresonant cap- 
ture processes, which is used in the NACRE compilation and many other studies, 
is a crude approximation, with Eq. flS}, of the accurate quantum three-body model 
calculation. The newly evaluated triple-a reaction rate will affect many studies on 
nuclear astrophysics, those on helium burning in accreting white dwarfs and neutron 
stars in particular. Detailed description of the theoretical framework together with 
further discussion on the comparison with other existing methods will be presented 
in a forthcoming paper. 

The authors wish to thank M. Kawai, M. Hashimoto, and E. Hiyama for helpful 
discussions. 
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